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ABSTRACT 

If cosmic ray protons interact with gas at roughly the mean density of the interstellar medium in starburst 
galaxies, then pion decay in starbursts is likely to contribute significantly to the diffuse extra-galactic back- 
ground in both 7-rays and high energy neutrinos. We describe the assumptions that lead to this conclusion and 
clarify the difference between our estimates and those of Stecker (2006). Detection of a single starburst by 
GLAST would confirm the significant contribution of starburst galaxies to the extra-galactic neutrino and 7-ray 
backgrounds. 

Subject headings: galaxies:starburst — gamma rays:theory, observations — cosmology:diffuse radiation — 
ISMxosmic rays — radiation mechanisms:non-fhermal 



1 . INTRODUCTION & PROTON CALORIMETRY 

Inelastic proton-proton collisions between cosmic rays and 
ambient nuclei produce charged and neutral pions, which sub- 
sequently decay into electrons, positrons, neutrinos, and 7- 
rays. Neutral pion decay accounts for most of the 7-ray emis- 
sion from the Milky Way at GeV energies (e.g., Strong et 
al. 2004). The contribution of pion decay to high energy emis- 
sion from other galaxies — and thus also to the extra-galactic 
neutrino and 7-ray backgrounds — remains uncertain, how- 
ever, because of the current lack of strong observational con- 
straints. For example, aside from the LMC (Sreekumar et 
al. 1992), EGRET did not detect 7-ray emission from star 
formation in any other galaxy. In two recent papers (Loeb 

6 Waxman 2006; Thompson, Quataert, & Waxman 2006; 
hereafter LW and TQW, respectively), we have estimated that 
starburst galaxies contribute significantly to the extra-galactic 
neutrino and 7-ray backgrounds, with most of the contribu- 
tion arising from starburst galaxies at z ~ 1 . TQW also pro- 
vide detailed predictions for the 7-ray fluxes from nearby star- 
bursts, several of which should be detectable with GLAST. 
Using what seem to be similar arguments, however, Stecker 
(2006) presents an upper limit on the starburst contribution 
to the neutrino and 7-ray backgrounds that is a factor of ~5 
below the predictions of LW and TQW. Here we clarify the 
difference between these predictions (see Sj3jl. 

The timescale for cosmic ray protons to lose energy via 
p-p collisions with gas of number density n cm" 3 is t pp w 

7 x 10 7 /nyr. The importance of pion losses depends on the 
ratio of t pp to the cosmic -ray proton escape timescale, r esc . 
The escape timescale for cosmic rays is quite uncertain, but 
in starbursts it is probably set by advection in a galactic su- 
perwind, rather than diffusion, as in the Galaxy. If correct, 
the escape timescale from a starburst with a wind of veloc- 
ity V « 300 km s" 1 and a gas scale height of h as 100 pc is 
of order r esc w h/V ~ 3 x 10 5 yr. This implies that if cos- 
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mic rays interact with gas of mean density n > 100 cm" 3 , then 
Tpp < T esc and we expect that essentially all cosmic ray pro- 
tons are converted into charged and neutral pions and their 
secondaries before escaping. In this limit, the galaxy is said 
to be a "cosmic ray proton calorimeter." If escape is instead 
due to diffusion rather than advection, then the density thresh- 
old required for a galaxy to be a proton calorimeter is likely 
to be significantly lower. 

In the Milky Way, few GeV cosmic -ray protons only lose 
k 10% of their energy to pion production before escaping 
the galaxy. For this reason, LW and TQW focused on the 
high energy emission from luminous starbursts. These galax- 
ies have much higher gas densities and likely dominate the 
contribution of star-forming galaxies to the total 7-ray and 
neutrino backgrounds. For reference, we note that the local 
starbursts NGC 253 and M 82 have n > 400 cm" 3 , and the nu- 
clei of Arp220 have gas densities of order 10 4 cm" 3 , well in 
excess of that required by the above arguments for t pp < r esc . 

Since the fraction of the proton energy supplied to sec- 
ondary electrons, positrons, neutrinos and 7-rays is deter- 
mined by the microphysics of the Standard Model, the 7-ray 
and neutrino fluxes from starbursts depend primarily on astro- 
physical properties of star forming galaxies via 

1. the fraction of a SN's canonical 10 51 ergs of energy sup- 
plied to cosmic ray protons, r\, 

2. the spectral index p of the injected cosmic ray proton 
distribution, 

3. the evolution of the comoving star formation rate den- 
sity of the universe, p\(z) (M Q yr" 1 Mpc" 3 ), 

4. and the fraction of all star formation that occurs in the 
proton calorimeter limit as a function of cosmic time, 
f(z). 

For any r\ and p, the 7-ray and neutrino emission from a 
given proton calorimeter can be readily calculated, while the 
cumulative extra-galactic backgrounds also depend on the star 
formation history of the universe via p*(z) and f(z). 

As TQW and LW show, in the proton calorimeter limit the 
total 7-ray and neutrino luminosities for an individual star- 
burst can be related to Ltir, the galaxy's total IR luminosity, 
by 

L 7 « (2/3)L nu w 1.5 X 10"Vo5£tir, (1) 
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FIG. 1. — The comoving star formation rate density p*(z) for a num- 
ber of models and fits to existing data: RSF1 (dashed black), RSF2 (solid 
black), and RSF3 (dotted black) are from Porciani & Madau (2001); the "Fos- 
sil" model (solid red line) is a simple broken power-law fit to Nagamine et 
al. (2006) (their Fig. 5); the "GALEX" model (solid blue line) is the dust- 
corrected p±(z), adjusted to the same assumed IMF as the other models pre- 
sented here, from Figure 5 of Schiminovich et al. (2005) (p*(z) assumed con- 
stant for z > 3); the "H&S" model (solid green line) is an analytic func- 
tion from Hernquist & Springel (2003). The total TIR background for these 
models as computed using eq. |2j is F TIR = 19.7, 20.4, 23.0, 15.9, 19.5, and 
17.3 nW nr 2 sT 1 for RSF1, RSF2, RSF3, Fossil, GALEX, and H&S models, 
respectively. These numbers for F TIR can be compared with observational 
constraints that indicate F TIR ~ 24- 27.5 nW mT 1 sr 1 (Dole et al. 2006). 
Taking f(z) = min(0.9z + 0.1, 1), / ca i = 0.79, 0.81, 0.81 0.76, 0.75, and 0.73 
for these same models. Using these values of f T1R and / ca |, the total 7-ray 
and neutrino backgrounds can be calculated from eq. Q as in eq. |4j. Figure 
|2| shows the cumulative fraction of F TIR as a function of redshift, £(z), for 
each of these models (eq.l3*l. 



where 770.05 = 77/O.O5 (i.e., 5 x 10 49 770.05 ergs per supernova 
is supplied to cosmic ray protons). ?/ is constrained by the 
requirement that the radio emission from secondary elec- 
trons/positrons in starbursts be consistent with the observed 
FIR-radio correlation (see, e.g., Yun et al. 2001). For 
proton calorimeters, 7/ > 0.05 is only possible if ioniza- 
tion, bremsstrahlung, and/or inverse-Compton losses domi- 
nate synchrotron losses for secondary electrons and positrons 
in starbursts (see TQW; Thompson et al. 2006). 

2. THE STARBURST BACKGROUND 

Because of the predicted one-to-one correspondence be- 
tween the neutrino, 7-ray, and radio emission of starbursts and 
their IR emission, the diffuse background in 7-rays or neutri- 
nos can be normalized to the TIR background, F T1R . Figure 
^ shows a number of models for p\(z) drawn from the litera- 
ture, based on observations of star-forming galaxies at differ- 
ent redshifts. The integrated background from star formation 
is then simply 



f tir = 



ep\c 2 dz 
4ttH J (1+z) 2 £Ti) : 



FIG. 2. — The cumulative contribution to the extra-galactic background 
in starlight ((z) (see eq.|3j as a function of redshift for the models of p*(z) 
shown in FigureHI For reference, the horizontal magenta lines have been put 
at £ = 0.4 and C = 0.5, and the vertical line at z, = 1. Note that in all models 
roughly half of the extra-galactic background comes from z > 1 ■ 



where E(z) = [ft m (l +z) 3 + ft A ] 1/2 and e is an IMF-dependent 
constant. Results for F T1R for the various models shown in 
Figure ^ are given in the caption and range from F™ « 0.8 
to F™ « 1.15, where F™ = F T1R /2Q nW m" 2 sr" 1 . Note that 
our calculation here does not include the order unity contribu- 
tion to the optical/NIR background at z ~ from the old stellar 
population (e.g., Nagamine et al. 2006; Dole et al. 2006). 

Figure [2] shows the contribution to F T1R as a function of 
redshift: 



C(z) = 



P\(z')dz' 



tJ 



P\(z')dz' 



(3) 



J (1+zOW)' [Jo (l + z') 2 E(z>) 
The dotted lines indicate z = 1 and £ = 0.4 and 0.5, and have 
been added for reference. Note that all of the models consis- 
tent with the observed star formation history of the universe 
have ~ 1 /2 of the TIR background produced at z > 1 . 

The total 7-ray and neutrino backgrounds can now be sim- 
ply estimated from the TIR background. Equation Q for the 
7-ray and neutrino emission from a given proton calorimeter 
implies that 

F 7 « (2/3)F na a 1.4 x 10^.05/0.75^ GeV s" 1 cm" 2 sr" 1 

(4) 

where / . 75 = / ca i/0.75 and 



p\(z)f(z)dz 



P*® dz \ (5 ) 
[(l+z) 2 E(z)]j 



(2) 



,o [(l+z) 2 £(z)]' uo 
is the fraction of the TIR background produced by starburst 
galaxies that are in the proton calorimeter limit (i.e., t pp < 
7~esc; ffl and f(z) is that fraction at every z. For a flat 7-ray 
spectrum (p « 2), equation © implies a specific intensity at 
GeV energies of 

i/i^GeV)* (2/3)1//^ 

w lO^WbsF™ GeVs-'cm^sr" 1 . (6) 



3 



Given the strong constraint on r/ implied by the FIR-radio 
correlation and the reasonable convergence of models and 
data on the TIR background (f 2 o IR w Fig-[D> the estimated 
7-ray and neutrino backgrounds depend primarily on the frac- 
tion of star formation that occurs in the proton calorimeter 
limit: f(z) and / ca i. Locally, we estimate f(z « 0) ~ 0.1 
from the fraction of the local FIR and radio luminosity den- 
sity produced by starbursts (e.g., Yun et al. 2001). However, 
at high redshift, a much larger fraction of star formation oc- 
curs in high surface density systems that are likely to be pro- 
ton calorimeters. For example, the strong evolution of the IR 
luminosity function with redshift implies that f(z) increases 
dramatically from z ~ to z ~ 1 (e.g., Dole et al. 2006). If we 
take f(z > 1) ~ 1 and pjg) as in any of the models of Figure 
^ we find that f ca \ ~ 1, with a significant contribution to the 
7-ray and neutrino backgrounds coming from z ~ 1-2. In 
fact, we find that 

F 1 oc i// tf (GeV) oc f(z « 1 - 2). (7) 

More specifically, taking Q, m = 0.3, f^A = 0.7, and assum- 
ing a function f(z) = min(0.9z+0.1, 1) that smoothly inter- 
polates from a small local starburst fraction (f(z s=s 0) s=s 0.1) 
to an order unity starburst fraction at high redshift, we find 
that / C ai ~ 0.8 for the models shown in Figure ^ ( see cap- 
tion). As an example of another alternative, assuming f(z) = 
min[0.1(l +z) 3 ,0.8], in which the starburst fraction increases 
roughly in proportion to p\(z) and saturates at f(z = 1) = 0.8, 
/cai ~ 0.6 for the same models. 

For numerical calculations for the 7-ray background spe- 
cific to a given f(z), p\(z), and proton injection spectra, p, see 
Figure 1 of TQW. 

3. COMPARISON WITH STECKER (2006) 

In a recent paper, Stecker (2006) (S06vl, v2, v3) estimates 
a significantly smaller 7-ray and neutrino background from 
starburst galaxies than found in TQW and LW, respectively. 
In his original version of the paper (S06vl), Stecker normal- 
ized the total neutrino background to the local FIR luminosity 
density from starburst galaxies in Yun et al. (2001) (based on 
IRAS), thereby neglecting the contribution to the background 
from redshifts above z ~ 0. This error led to a factor of wlO 
underestimate of the total neutrino background relative to LW. 

S06v2 corrects this error and addresses both the 7-ray 
and neutrino backgrounds expected from starburst galaxies. 
Stecker finds a value for the neutrino and 7-ray backgrounds 
approximately equal to the value quoted in S06vl, ^5 times 
lower than LW and TQW. There are two primary differences 
that lead to this discrepancy: 

1. Stecker advocates that 22% of the extra-galactic back- 
ground is produced by star formation above z > 1.2 (see 
his Table 1, S06v3). Figures ^ an d ID show, however, that 
essentially all models for the star formation history of the 
universe that are consistent with the observed evolution of the 



star formation rate density with redshift (p\(z)) yield roughly 
equal contributions to F T1R below and above jfsl, 

2. We argue that order unity of all star formation at z ~ 1 
occurs in starbursts likely to be proton calorimeters, whereas 
Stecker quotes 13% in the range 0.2 < z < 1.2 and 60% for 
z > 1.2 (see his Table 1, S06v3). This step function model 
for f(z) likely significantly underestimates the contribution 
of z ~ 1 galaxies to the 7-ray and neutrino backgrounds. 

Thus, the inconsistency between our predicted 7-ray and 
neutrino backgrounds and Stecker's upper limit is a simple 
consequence of his model for the star formation history of the 
universe and the redshift evolution of the starburst population. 
For the reasons explained above, we believe that our assump- 
tions are in better agreement with observations of star-forming 
galaxies at z > 1 . 

4. CONCLUSION 

The magnitude of the observed diffuse 7-ray background 
is quite uncertain, primarily as a result of foreground 
subtraction. Whereas Sreekumar et al. (1998) (compare 
with Strong et al. 2004) find a specific intensity of 1 .4 x 
10" 6 GeV s" 1 cm" 2 sr" 1 at GeV energies, Keshet et al. (2004) 
derive an upper limit of vl v < 5 x 10~ 7 GeV s" 1 cm" 2 sr" 1 . 
Future observations by GLAST should alleviate some of the 
uncertainty in the background determination. A comparison 
of the upper limit by Keshet et al. (2004) and equation l|6) 
indicates that starburst galaxies may indeed contribute signif- 
icantly to the diffuse 7-ray background. 

Even in stacking searches, EGRET did not detect a single 
starburst galaxy (Cillis et al. 2005). In TQW (their Figure 1) 
we show that the EGRET non-detections are required by pre- 
dictions that are correctly calibrated to the FIR or radio emis- 
sion from galaxies. Hence, the EGRET non-detections do not 
provide a significant constraint on the contribution of starburst 
galaxies to the 7-ray or neutrino backgrounds. GLAST will, 
however, be much more sensitive. Using equation Q and nu- 
merical calculations for different p, TQW show that GLAST 
should be able to detect several nearby starbursts (NGC 253, 
M82, IC342; see TQW Table 1) if they are indeed proton 
calorimeters. 

The dominant uncertainty in assessing the starburst contri- 
bution to the 7-ray and neutrino backgrounds lies in whether a 
significant fraction of cosmic -ray protons in fact interact with 
gas of average density or whether the cosmic rays escape in a 
galactic wind without coupling to the high density ISM (the 
star formation history of the universe is comparatively well 
understood). Detection of a single starburst galaxy by GLAST 
would strongly support the estimates of LW and TQW and 
thus the importance of starburst galaxies for the extra-galactic 
neutrino and 7-ray backgrounds. Such a detection would also 
provide an important constraint on the physics of the ISM in 
starburst galaxies. 
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